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Abstract 

We present theoretical results on microwave and far-infrared (FIR) absorption of single-electron 
transistors obtained within exact numerical diagonalization for finite clusters. They show that both 
the microwave and the FIR spectra consist of two maxima, whose origin can be understood physi- 
cally. Our results on microwave absorption provide a physically intuitive qualitative interpretation 
of the Kondo splitting observed by Kogan et al [Science 304, 1293 (2004)]. The present results 
on the FIR absorption supplement and provide a physical insight into previous results obtained 
by means of the numerical renormalization group. Based on our theoretical results, we propose to 
conduct FIR experiments to determine the charging energy and other relevant parameters. 

PACS numbers: 85.35.Gv, 36.20.Ng, 73.23.Hk, 81.07.Ta, 73.63.Kv 
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I. INTRODUCTION 



In a single-electron transistor (SET), which consists of a quantum dot (QD) attached to 
two electrodes, a small source-drain voltage yields a current flowing only for certain values 
of the gate potential Vg.-'^^ At temperatures below the Kondo temperature (T < T^), 
conduction occurs in a V^-range delimited by the situations where the energy of the dot 
level Ed is such that the lowest Hubbard "band" or the highest Hubbard "band" are nearly 
resonant with the electrode Fermi energy ep, Sd ~ and + f/ ~ Sp, respectively. Here, 
U represents the dot charging energy, i. e., the energy required to add an extra electron on 
the dot. The zero-bias conductance G reaches the unitary limit G = Gq = 2e^ /h within the 
Kondo plateau, which is deflned hy ep — U ^Sd^ sp. 

The charging energy represents a key parameter for SETs. The unpleasant fact is that in 
dc transport measurements of the zero-bias conductance U cannot be directly determined, 
because the dot energy Sd cannot be directly controlled, but rather indirectly via the potential 
of a "plunger" gate potential Vg, on which it linearly depends: Sd = aVg + const. If Vg^i and 
Vg^u denote the gate potentials whereat the lower and upper Hubbard bands become resonant, 
one gets U = aiVg^i — Vg^)- So, to determine U, in addition to the difference Vgj — Vg^u, 
which is available from the transport data, supplementary hypotheses are needed to deduce 
the conversion factor a, which, although physically plausible, cannot be fully justifled at the 
nanoscale and require assumptions or arguable extrapolations of macroscopic relations to the 
nanoscale. One way is to assume a certain phenomenological T-dependence (convolution of 
a Lorentzian with the derivative of the Fermi function) and to flt the width of the Coulomb 
blockade peaks G{T)^ Another possibility is to resort to the capacitance model, which 
describes the SET in terms of three effective capacities Gg, Gg, and Gd between the dot and 
the gate, source, and drain, respectively.-i^>^ The conversion factor, expressed by a = Gg/G 
[G = Gg + Gg + Gd), can then be obtained from the Coulomb diamonds of the stability 
diagram. Even without inquiring whether such assumptions are justifled, the inaccuracies 
of the parameters estimated in this way are rather large; uncertainties can be as large as 
~ 20%.^ Therefore, utilizing more accurate or at least alternative methods of investigation 
is highly desirable. 

It is a main goal of this paper to show that the far-infrared (FIR) absorption represents 
a possible alternative technique for the characterization of SETs and how such experiments 
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can be conducted. 

The remaining part of this paper is organized in the following manner. In Sect. [TTl we 
expose the theoretical framework, and in Sect. IIIII present all relevant computational details. 
In Sect. IIV[ exact numerical results for the full ac absorption spectra are presented and 
analyzed in terms of a few physically relevant many-electron states. Sect. |V] is devoted 
to finite size effects. Next we discuss the two distinct spectral ranges significant for SETs 
separately: the microwave/radiofrequency absorption in Sect. |VT]and the FIR absorption in 
Sect. IVIII Experimental implications of the theoretical results for the FIR absorption are 
presented in Sect. IVIIII Sect. IIX]is devoted to conclusions. 



II. THEORETICAL FRAMEWORK 



Following the usual procedure, we shall describe the SET within the Anderson single- 
impurity model^'^'^'^ 
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The left (L) and right (R) electrodes are assumed to contain noninteracting electrons, which 
are characterized by the same bandwidth D = At and the same coupling td to the dot. The 
dot is modeled by a single level, whose energy Ed can be tuned by means of a gate potential, as 
discussed in Sect. HI a„_o- (dn^) are annihilation (creation) operators for electrons in the left 
and right leads {L,R) and d^j = a^^^ (dj. = aj^) destroys (creates) electrons in the QD. The 
number of electrons will be assumed to be equal to the number of sites, N = Ml + Mn + 1. 

The quantity of interest, the frequency-dependent absorption coefficient /i(co') in the 
ground state \l/o (case of zero temperature), can be expressed as a sum of contributions 
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of various excited states \I/a (-f^^A = -E'a^E'a 



(2) 



where Pd is the QD dipole moment. Eq. ([2]) represents the result of the hnear response theory 
by considering an ac electromagnetic perturbation H' = —P^So cos out. Various aspects of 
the problem of a SET in an ac field within the linear response approximation were previously 
considered in several studies (see, e. g., Refs.QQQ. Because the definition of the dipole 
operator P^ for a point-like QD poses some problems, it is more convenient to express the 
matrix elements entering Eq. ([2]) (\l'A|-Pd|^o) = — '^^(^A|jd|^o)/(-E'A — Eq) in terms of the 
current operator jd = (j_i/2+Ji/2)/2, as done in similar cases,~ which can be unambiguously 
defined as jn+1/2 = itn{e/h) J2a(^n a^n+i,a — H.c). So, the ac absorption is specified by the 
spectral lines A characterized by the absorption intensities fi\ and the absorption frequencies 
ujx defined by 
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fix = — |(^A|r|^o) 

UJx 



(3) 



In the presentation and the discussion of the results on ac absorption, unless otherwise 
specified, we shall refer throughout to a SET in the Kondo regime {sp — U ^ £d ^ £f)- 
Moreover, we can restrict ourselves to the range Sp ~U /2 < Sd sp because of the particle- 
hole symmetry. 

III. COMPUTATIONAL DETAILS 



Below, we shall present results on the ac absorption of a SET obtained by exact (Lanczos) 
numerical diagonalization. The method of com put ation we employ here is that used in our 



earlier works; see, e. g., Refs. 
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Because the full details on this 



method were not published and because of the significant differences between our Lanczos 
implementation to compute the linear response and the more familiar continued fraction 
algorithm,— we describe it below for the benefit of the reader. 

In the first run, the Lanczos procedure is iterated until, after Mc iterations, the lowest 
(ground state) energy Eq converges. In the second run, by carrying out again J\fc iterations 



and with the same starting Lanczos vector, the corresponding Ritz vector \I/o is computed 
by accumulation without the need of storing the Lanczos vectors. To check that this vector 
represents indeed the accurately evaluated ground state \I'o, we straightforwardly compute 
the dispersion (\l'o|(// — -E'o)^|^o)^^^ and convince ourselves that it is much smaller (usually 
5 — 6 orders of magnitude) than the lowest excitation energy. The above scheme can also 
be used to reliably compute several lower excited eigenstates, but it is usually unpractical 
to target all the eigenstates \E'a needed to compute the linear response via Eq. ([2]), e. g., by 
orthogonalization on eigenvectors already converged in previous runs. The reason is that 
many eigenvectors, which are not important for the linear response, are also targeted. To 
ensure that the important eigenvectors are targeted, in a third Lanczos run, we employ a 
starting Lanczos vector adequate for the specific linear response considered. This is, in the 
present case, the normalized vector Pdl'^o)- The needed matrix elements (\E'A|-Pd|^o) are 
given by the first component of the tridiagonal vectors obtained in this third run. Usually, 
a number of iterations comparable to JVc suffices for the third run. As an important test 
of the results for the linear response computed in this way, we always check whether they 
satisfy the sum rule, which can be deduced exactly from Eq. ([2]) 

J]|(vl/,|P,|vl/o)|2 = (vl>o|p]|^o) , (4) 

A 

because the r.h.s. is known, namely the squared norm of the vector Prf|\E'o). In certain cases, 
the linear response computed within the third run does not satisfy the above sum rule, e. g., 
because of spurious vector duplication. Therefore, to be always on the safe side, we carry out 
a fourth Lanczos run, wherein, similar to the second run, we also compute and store all those 
Ritz vectors \E'a, which where found to have a significant spectral weight | (\E'A|Pd|\E'o) P [in 
practice, above 10~^ of the r.h.s. of (Eq. (jlj)] in the third run. Storing these vectors \E'a is not 
much more demanding than storing the ground state \E'o alone, because for all the problems 
we investigated so far, at most ~ 10 — 20 Ritz vectors are important. The real, prohibitive 
limitation remains, as in all exact diagonalization approaches, the cluster size. Again, we 
check that these Ritz vectors are accurate eigenvectors by straightforwardly computing the 
dispersions {^x\iH — _E'a)^|^a)^^^- By using these eigenvectors \E'a we finally compute the 
linear response from Eq. ([2]) and convince ourselves that all important eigenvectors have 
been targeted by checking the sum rule (jlj). 

Proceeding in this way, the computing time is at most ~ 1.5 — 2 times larger than for 
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implementations of the continuous fraction algorithm,— but we can safely rule out any 
numerical artefacts and have the guarantee that the solution obtained is mathematically 
exact. In addition and equally important, this method allows us to compute and resolve 
individual nearly degenerate spectral lines, a situation where the information that can be 
extracted from convoluted spectra provided by the continued fraction algorithm does not 
suffice. This represents a quite relevant aspect for SETs and other QD-based nanosystems, 
where nearly degenerate states with the same symmetry (avoided crossings) are often en- 



countered; see Refs. 
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m and Sect. |VT1 



IV. EXACT RESULTS ON THE FULL AC ABSORPTION SPECTRA AND 
THEIR PHYSICAL INTERPRETATION 

Numerical exact results for frequencies and intensities of all the ac absorption signals 
obtained as described in Sect. IIIII are collected in Fig. [T] They have been obtained for 

= 11 and parameter values, which are typical for real cases: t = 0.5 eV (electrode 
bandwidth D = At = 2eV), = 0.2 meV, and U = 8meV. We emphasize that these are 
numerical exact results, obtained by using all the 213444 multielectronic configurations of 
the eleven-site cluster with eleven electrons and a total spin projection Sz = 1/2. Based on 
the considerations of Sect. IIIII we can safely state that the ac spectrum of the investigated 
cluster solely consists of four relevant absorption signals. The other transitions, although 
allowed by symmetry, are completely irrelevant, as their intensities are orders of magnitude 
smaller and are therefore invisible in Fig. [Dd.— 

Exact results on the SET ac absorption have of course their own importance, but do 
not yet provide much physical insight into the problem. Since the above exact results show 
that only four optical transitions are important, one can expect that, out of numerous 
multielectronic configurations (namely, 213444, see above, in the case under consideration, 
of an eleven-site cluster with a total spin projection Sz = 1/2), there should only exist a few 
many-body states, which are relevant. If so, the problem is of course to identify them and 
to unravel their physical content. This shall be done next. 

There are nine such significant many-electron states. These configurations (|1) to |9)) are 
schematically shown in Fig. [21 Configurations |1) to |5) correspond to one electron on the 
dot, in |6) and |7) the dot level is vacant, while in |8) and |9) it is occupied by two electrons. 
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FIG. 1: (Color online) Dependence on the dot energy of the four relevant (a) absorption fre- 
quencies ijj (in meV) and (b) absorption intensities [i (in arbitrary units) of the optical transitions 
with significant spectral intensities. In (a), deeper within the Kondo plateau, the exact frequencies 
a;3,4 are well approximated by The dashed lines are the analytical continuations of the full 

lines. The parameter values are: t = 0.5 eV, td = 0.2 meV and U = 8meV. 



A superficial glance at the schematic representation of Fig. |2] can easily overlook both the 
underlying physics and the computational effort involved, and therefore a comment is in 
order at this point. Out of the electrons in the two electrodes, only those occupying the 
Fermi levels are shown for the nine states of Fig. [2l For these nine states, the single-particle 
states of the electrons in the electrodes are in momentum (k) space, and not in the real (site, 
n) space, in which the exact numerical diagonalization is carried out because the Hamiltonian 
matrix, Eq. ([T]), is sparse. A single-particle fc-state, e. g., in the left electrode represents a 
superposition of single-particle n-states. In addition, one should note that the electrons 
in electrodes depicted in Fig. [2] represent electrons at the Fermi level. This means that these 
electrons are delocalized over the electrodes. Consequently, although we show below that 
the approximative description in terms of the nine relevant states is accurate, it is not a 
priori obvious that the problem can be reduced or reasonably approximated by studying 
a three-site cluster. To summarize, each of the nine states depicted in Fig. [2] contains in 
fact numerous multielectronic configurations in the real space. However, what is physically 
important is the existence of a very reduced number of the relevant states. 

The discussion below proceeds in terms of these nine many-body states with significant 
contributions to the ground state and the four excited states 1, 2, 3, and 4 depicted in Fig. 
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FIG. 2: (Color online) Multielectronic configurations with significant contributions to the ground 
state ^'o and the excited states ^'1,2,3.4 important for ac absorption. For each configuration, we 
show the electrons at the Fermi levels of the left and right electrodes, and on the dot (red, blue, 
and green, respectively). In either electrode, the single-electron states below the Fermi level are 
occupied. 

[TJ From these nine most relevant states one can construct the following states with definite 
spin Sz = S = 1/2 (notice that the total electron number is odd), which are either even (^f) 
or odd (m) under space inversion 
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The eigenstates important for ac absorption can be well approximated as 

l^o) - cosx - 1^2) sinx 

1^1) ^ \gi) cos e -Igs) sine ^ \gi); 

1^2) ^ 1^72); (6) 
1^3) ^ |5'i)sin6'+ \gs)cose ^ Ig^); 
1*4) ^ 1^74). 

Eqs. hold for ep — U/2 < Ed < Sp- We can restrict ourselves to this range because of 
the particle-hole symmetry. For sp — U < Sd < Sp — U the states |6) and |7) must be 
replaced by |9) and |8), and vice versa. 

To illustrate that the eigenstates ^1^0,1,2,3,4 computed exactly are indeed very well approx- 
imated by the expressions in the r.h.s. of the symbols ~ in Eqs. (El), we present in Figs. [3] 
and m the curves of the weights = l(^i,2|^o)P and = |((7j|\l'j)p = 1,3). In all 
these cases, the two functions entering the expressions in the r.h.s. of Eqs. ([6]) exhaust the 
expansions of the exact eigenstates \l'o,i.3 within an accuracy of ~ 10~^. This fact fully jus- 
tifies the use of the intuitive notations in terms of cosines and sines in Eqs. (|6]), cos^ x = Pi) 
cos^ 9 = p\. As concerns the other two exact eigenstates, the approximations |^t'2,4) — |fi'2,4)5 
are also accurate within ~ 10~^. As visible in Figs. [3] and HI deeper within the Kondo 
regime, cosx — 1 and cosO ~ 1, and therefore 1^^0,1,3) are reasonably approximated as ex- 
pressed in the r.h.s. of the arrows in Eqs. ([6]). Bearing this in mind and inspecting Eqs. ([6]) 
and ([5]) and Fig. [21 one can identify two groups of important eigenstates, which are well 
separated energetically. The first group comprises the eigenstates ^1^0,1,2; which basically 
consist of superpositions of the nearly degenerate configurations |1) — |5), corresponding to 
states with a singly occupied dot. This fact nicely reveals the spin entanglement and the role 
of the coherent superpositions of all the possible spin flip processes (|1) ^ |3), |3) ^ |5), 
|2) ^ |3), |4) ^ |3)) in the formation of the nearly degenerate states \t'o,i,2 important for 
the Kondo effect. The absorption frequencies uJi^2 of these optical transitions are low, falling 
into the microwave-^^ or even radiofrequency (rf) range. The second group comprises the 
higher energy states \t'3 and \t'4, which correspond to a dot that is either doubly occupied 
or empty. Loosely speaking, they amount to excite a particle-hole pair, wherein the hole 
state is on the dot and the particle state in electrodes, or vice versa. The corresponding 
absorption frequencies, uj^ ui^ = ep — Ed and a;4 ~ a;^ = + ^ — (cf- Fig. [1^), are of 
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FIG. 3: (Color online) e^-dependence of the weights 2 = K'Wi,2|'J'o)P of the states ui^2 entering 
the linear combination of Eq. ^ for clusters with = 3, 7, 11 sites. Parameter values as in Fig. [TJ 
For all A^'s, the deviation from unity of the sum p\ + P2 (a-t most ~ 10~^) is invisible within the 
drawing accuracy. 
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FIG. 4: (Color online) Curves for = |(5ri|vl'i)p and pg = 1(5(3 l^fi)]^ similar to Fig. [3l For all TV's, 
the p{-curve cannot be distinguished within the drawing accuracy from that for = |((73|^'3)p, 
and the pg-curve from that for pf = | ((jilt's) p. For all A^'s, the deviation from unity of the sum 
Pi +P3 (at most ~ 10^'^) is invisible within the drawing accuracy. 

the order of the charging energy U, falling therefore into the FIR range. 
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V. FINITE-SIZE EFFECTS 



As is well known, the drastic limitation of the exact numerical diagonalization to rather 
small cluster sizes N often precludes a reliable finite scaling analysis. There are well known 
examples (see, e. g., Refs. Il7|j32|j33l ) of non-monotonic A^-dependent properties, or qualita- 



tively different behaviors at smaller and larger due to a different underlying physics (see. 



e. g., Ref. 



33l ) at the sizes where exact numerical diagonalization is feasible. This limitation 



is even more severe in the case of SETs, in the sense that not even all these small sizes 
can be included in a finite-scale analysis. A careful selection of the A^-values to be included 
in the finite-scale analysis is often necessary, as is well known, e. g., in the case of cyclic 
polyenes CjvHtv or related systems, where Hiickel (A^ = An + 2) and anti-Hiickel (A^ = 4n) 
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34 and references cited 



systems behave differently (n is an integer); see, e. g., Refs. 19||i 
therein. With our implementation described in Sect. IIIIl we can reliably treat the linear 
response of half-filled clusters up to A^ = 14, amounting to a dimension of the Hilbert space 
of 11,778,624. This is not too much different from the largest size (A^ = 12) of most recent 
studies on the dc-conductivity of model (^.^^ In view of the analysis in terms of the relevant 
many-body states of Fig. [21 it is clear that considering symmetric clusters (identical elec- 
trodes) is advantageous. Because short electrodes with an even number of sites are known to 
yield spurious results (compare Ref. 
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with Ref. l35l). what remains is to consider electrodes 
with an odd number of sites, which mimic "metallic" electrodes (i. e., electrodes with a 
partially occupied Fermi level).— Concretely, this means that we are left with the values 
A^ = 3, 7, 11. Obviously, one cannot expect to rehably deduce a scaling law solely based on 
these three A^-values. 

In view of the aforementioned limitations, similar to our previous works,— »^ we shall 
simply inspect whether the relevant properties computed for A^ = 3, 7, 11 are significantly 
size dependent or not. Typical results are shown in Figs. [3], HI andO They reveal that certain 
quantities, like the lowest excitation energies u!i,2 of Fig. [5b are strongly size dependent. 
Obviously, such results for ui^2 of exact diagonalization cannot be trusted, at least not 
quantitatively (see also Sect. [VTl) . 



37 and 
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there also exist quantities. 



But, similarly to the examples presented in Refs. 
which only slightly depend on A^. Most important for the main purpose of this work, this 
is the case of the two higher absorption frequencies cus^ of Fig. [5K, the key quantities to 
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FIG. 5: (Color online) Results on the (a) higher (FIR) frequency a;3^4 and (b) lower (rf/microwave) 
frequency ui^2 ac absorption for several cluster sizes N and same parameter values as in Fig. [TJ In 
panel (a), the triangles and circles are for clusters with = 3 and N = 7, respectively, and solid 
lines are for clusters with = 11. The latter cannot be distinguished within the drawing accuracy 
from those of asymmetric clusters, wherein the dot is attached to a single electrode with 7, 9, 11, 
and 13 sites. In the Kondo regime, 003^4^ are only slightly size dependent, while u!i^2 are strongly 
size dependent. Notice the logarithmic scale on the ordinate in panel (b). 

be measured in the FIR experiments we propose here (see Sect. IVIIII) . Therefore, to give 
further support to the fact why we believe that, deeper in the Kondo plateau, the results 
for the curves uj3^4{ed) are not significantly affected by finite size effects, we carried out 
supplementary calculations. Namely, we considered asymmetric clusters, wherein the QD is 
attached to the end of a single "metallic" electrode with an odd number of sites A^„. This 
procedure, which amounts to unfold the original symmetric cluster,— i^i^^ has the advantage 
that the size of the single electrode can be larger, roughly twice that of one electrode of a 
symmetric cluster. The largest relevant (odd) size that we can treat by exact diagonalization 
is Nu = 13 (13 + 1 sites). The shortcoming of the asymmetric cluster is that it misses the two 
lowest excitations u)i^2 related to the coherent spin fluctuations responsible for the Kondo 
effect. The first excitation of the asymmetric cluster, which is almost degenerate with 
the singlet ground state, is a spin triplet, and the small singlet-triplet splitting could be 
considered as the counterpart of co'1.2 in symmetric clusters. However, this triplet excited 
state is irrelevant for the spin conserving ac absorption processes. Most important is that 
the next two excitations of the asymmetric cluster are singlet states, which are optically 



12 



active, and their energies are the counterpart of the above c<j3 4. As noted in the caption of 
Fig. [5^, the curves for Nu = 7, 9, 11, 13 cannot be distinguished from those of the symmetric 
cluster with = 11, which is the counterpart of the asymmetric cluster with Ny_ = 5. 

For completeness, we mention that the size dependence of 003^4 remains weak even beyond 
the Kondo plateau (cf. Fig. E]), although this fact is not very important because of the small 
absorption intensities (cf. Fig. Wp)- There, the physical character of the Ci;3^4-excitations 
is different. Within the Kondo plateau they are related to excitations of a particle-hole 
pair, while beyond the mixed valence points they are related to excitations of two particle- 
hole pairs. This becomes clear if one inspects Fig. [5^, where the energies of the latter 
processes in the absence of electrode-dot coupling {td 0) are represented by the thin lines 
2{ed -eF)+U = {Ed - sf) + (e^ + f/ - Ef) and 'l{ep -ei)-\J = {ep - Ed) + {sp -Ed- U)M 
A similar change in the physical character can be seen, e. g., in the mixed valence region 
between the singly occupied and the vacant dot. There, the curve u^, which corresponds 
to the excitation of an electron from the singly occupied dot into electrodes {uj^ ~ tUg = 
Ep — £d), evolves into that amounting to bring an electron from electrodes to the vacant dot 
{u3 ^ Ed — £f)', see the lower right corner of Fig. 

By inspecting Figs.[3]and|4l one may argue that the size dependence of the wave functions 
\['o, ^1, and \E'3 is comparable; so, where does the difference between the size dependence of 
uJi^2 on one side and ti;3,4 on the other side come from? The reason is the following. While 
the cji 2-values are close to 0, the to'3^4-values vary close to the 4-values. which correspond 
to electrode-dot excitations in the limit of vanishing electrode-dot coupling (td 0), and are 
large (~ U) deeper within the Kondo plateau. In fact, the size dependence of the difference 
Us^/i — 1^34 is comparable to that of u!i^2 — 0, as seen in Fig. [5^1. It is the same strong A^- 
dependence of uji^2 that makes /ii,2 [cf. Eq. (|3])] strongly size dependent; the matrix elements 
of the hopping operator f are nearly A^- independent. 

Although the results presented above have shown that the size dependence of the two 
higher optical transitions is not substantial, the important question is whether the absorption 
peaks /i3,4 survive when the cluster is linked to infinite electrodes. Based on our previous 
investigation of photoionization^^ and on extensive calculations of the FIR absorption in 
broad ranges of SET parameters, we expect the following. As the size increases, the single- 
electron levels in electrodes become more and more dense, and the straight lines of 
Figs. [U or will intersect the numerous horizontal (i. e., e^-independent) lines corresponding 
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to excitations of particle-hole pairs in electrodes, in a way similar to that of the energies of 
various ionization processes (one-hole, two-hole-one-particle, etc) shown in Fig. 2a of Ref. 



Similar to Ref. 



24 this gives rise to a sequence of avoided crossings, but the spectral intensity 



remains concentrated in two diabatic states, as if these intersections were absent. From this 
perspective, one can also understand why, sufficiently away from the mixed-valence points, 
the curves for u!3^4,{ed) for the symmetric three-site cluster of Fig. [5^ represent reasonable 
approximations: roughly, they correspond to one electron-hole pair excitations, wherein the 
state of one mate of the pair is on the dot and the other at the electrode Fermi level. 
For the same reason, even the asymmetric two-site cluster provides a qualitatively correct 
description of the FIR absorption. 

As is well known, a weak electrode-dot coupling yields a small broadening (F ~ 2t^/t) 
of the isolated dot level Ed- The analysis of Sect. [IV] indicated that, basically, each of these 
transitions amounts to excite an electron-hole pair. Therefore, the electrode-dot coupling 
should reflect itself in a small broadening of the FIR peaks centered on the values — 
uj^ = Ef — Sd and oj^ ~ uj^ = Ed + U — ep, which replace the delta-shaped /xs^-lines of 
the finite cluster. From a strictly mathematical standpoint, to demonstrate that these FIR 
peaks survive when the finite cluster is linked to real electrodes, we can simply invoke 
their presence in the numerical renormalization group (NRG) results,-^ which are exact and 
consider infinite electrodes. 



VI. RADIOFREQUENCY/MICROWAVE ABSORPTION 

The existence of two electromagnetic transitions |\E'o) I ^1,2) with low absorption fre- 
quencies uJi^2 in the rf/microwave range is a remarkable theoretical result, because it is 
directly related to the recent experimental findings in SETs irradiated with microwaves.— 1^ 
Unfortunately, at present we cannot offer a reliable quantitative analysis and must restrict 
ourselves to a few qualitative considerations. The first, obvious reason of this impossibility is 
the strong size dependence of the results discussed in Sect. |Vl But there still exists another 
reason. As the electrodes become longer and longer (A^ 00), we expect that ui tends 
to the width of the Kondo resonance ~ Tk- At larger U, this width falls off exponentially 
with U, while our exact diagonalization data exhibit a much weaker, power law decrease 
with U. This [/-dependence is similar to that of the width in the density of states obtained 
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within a one-particle Green function approach.-^^ In that approach, also adopted in a series 



of other works (see Ref. 



35 



and citations therein), the finite cluster is embedded into infinite 
electrodes via a Dyson equation, wherein the self-energy is supposed to be not affected by 
electron correlations. We are not aware of similar developments for the two-particle Green 
function needed to compute the ac absorption. Still, the aforementioned similar and (in 
this respect) incorrect ^/-dependence of that approach and the present one seems to signal 
the need for a method that (presumably approximately but accurately enough) accounts for 
correlations in clusters of sizes much larger than the exact diagonalization can handle. In 
this sense, we think that the description of Sect. IIVI in terms of a few relevant many-body 
states is useful, since it emphasizes the similarity of the lowest two frequencies 001^2 to a tun- 
nel splitting. The coherent spin fluctuations embodied into the functions \l/o,i,2 expressed 
by Eqs. ([5l [6]) amount to a coherent tunneling between configurations that are classically 
degenerate and have indeed similarities to the tunneling between the degenerate minima of a 
symmetric double well potential. Most relevant, exponential decays of the tunnel splittings 
with the interaction strength are typical.—'^ In view of the severe size limitation within 
exact numerical diagonalization, and because it is unlikely that the small difference between 
Ui and UJ2, which becomes much smaller at larger sizes, can be resolved within the density 
matrix renormalization group (DMRG), we believe that at present the only possible ap- 
proach is a semi-analytical one, e. g., based upon symmetry-adapted trial wave functions for 
the lowest states \l/o,i,2, which also turned out useful for other strongly correlated electron 
systems.— 

To end this section, we believe, in spite of the above somewhat speculative considerations, 
that one can plausibly ascribe the excitation energy ui as the width of the Kondo resonance, 
while the excitation energy uj2, close to but still different from uJi, can be interpreted as the 
splitting of the Kondo resonance observed experimentally.— 

VII. FIR ABSORPTION 

In this section we shall focus on the other two transitions |\l/o) —>■ 1^3,4)- As seen in Fig. 
[IK, the absorption frequencies ^^3^4 are of the order of U. For many fabricated SETs (see. 



e. g., Refs. , and 



43I ) these values belong to the FIR range. The explicit forms and 



show that in the Kondo regime these two transitions amount to excite the electron from 
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the QD lower Hubbard band into the electrode Fermi level, and from the electrode Fermi 
level into the QD upper Hubbard band; sufficiently away from the mixed valence ranges 
{ed^epj^d^^F — U), the exact excitation energies are well approximated hj = ep — Ed 
and u'^ = ed + U — Ep (see Fig. [T^). 

Based on Fig. [H one expects in general two absorption peaks of a SET irradiated with 
FIR radiation. In the middle of the Kondo plateau {e^ = ep — U/2) the two transitions 3 
and 4 are degenerate, and therefore a single peak can be observed experimentally. There, 
the absorption frequency is just one half of the charging energy, = = u* ~ U/2. 
By moving away from this point in either direction, the absorption peak splits into two 
peaks of different intensities located symmetrically with respect to the degenerate peak, 
UJ3^4 ~ cij* =p \ed — The farther from the symmetric point, the more pronounced is the 
asymmetry in intensity, the stronger is the peak fi^ at the lower frequency 003 , and the weaker 
the peak at the higher frequency u^. 



Out of the studies on SETs in ac fields,— i^*^ excepting in part for Ref. 



12, none consid- 



ered the above aspects. Without establishing any relationship to the FIR absorption, the 



numerical results on frequency-dependent conductance deduced within the NRG of Ref. 



12 



show, interestingly, a weak peak (to which the authors paid little attention) for two values of 



12, to which 



Ed'- aX Ed = Sp — U and aX Ed = £p — U /2 (see Figs. 2 and 3, respectively of Ref. 
we refer below). This peak is directly related to our results. The situation Ed = Ep — U /2 
corresponds just to the point of particle-hole symmetry, and the peak position is visible, 
just as predicted by the present approach, at a; = a;* (note that Ep \s set to zero in Ref. Il2l). 



For Ed = Ep — U, the peak in Fig. 2 of Ref. 



12 



occurs at ~ (6 X IO-7O.O25) U = 0.2AU 



but the authors provide no physical interpretation of this value. In excellent agreement with 
this value, the lower frequency absorption peak /is predicted by our approach is = U/4. 
In addition, we predict another absorption peak /i4 at a higher frequency ^4 = {3/4)11 



UJ4/D = 0.01875 in the notation of Ref. Il2l). which, although in the range showed in 



'ig- 3, 
ia,ii we 



is invisible there. We can explain this fact: for the parameters employed in Ref. 
estimate that the higher frequency peak would be one order of magnitude less intense than 
the lower frequency one. This weak intensity could hardly be distinguished in the back- 
ground of the curve of Fig. 3 a.t u/D = 0.01875. To reveal the two peaks in FIR absorption, 
the NRG calculations should have used situations sufficiently away from the particle-hole 
symmetry point {\Ed — Ed\ should exceed the peak widths) but still sufficiently close to it. 
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because otherwise the high frequency peak would be too weak and thence not visible. 

To end this section, we note that the two peaks in the FIR absorption at uj-^ and UJ4, 
are the counterparts of two maxima located close to the energies Sd and Sd + t/, which are 



present in the electronic density of states (DOS~ 
to the Kondo resonance (see, e.g., Fig. 3 of Ref. 



ong with the sharp peak corresponding 



VIII. EXPERIMENTAL IMPLICATIONS 

Based on the above theoretical results, we propose to employ the FIR absorption as an 
experimental tool to characterize SETs. To avoid misunderstandings, we emphasize that 
the proposed experiments are different both from those carried out using rf or microwave 



radiation suitable for studying the Kondo resonance (e. g, Ref. l3ll ) and from those recently 
proposed by us to use photoionization,— where photons should have energies of the order of 
the work function (ultraviolet radiation). 

In experiments, even using a very well focused flux of FIR photons to irradiate a SET, it 
is important but, fortunately, easy to discriminate between absorption processes occurring in 
the dot, and in electrodes or due to acoustic phonons. One should simply monitor absorption 
by varying V^: the former signals are affected and should be analyzed, while the latter are 
not and should be disregarded. To exploit the present results, most desirable would be to 
record FIR absorption spectra of SETs directly. The absorption intensities may be very 
weak and their measurement a challenge for experimentalists. Even though difficult, this 
can no longer be considered a hopeless experimental task, particularly in view of the very 
recent advances in the field of molecular devices, enabling to measure the photon emission^ 
or Raman response^^ of a single molecule. 

As an easier experimental task, similarly to our earlier work^^, we propose to perform 
a mixed FIR-absorption-dc-transport study, which should not pose special experimental 
problems. Again, the fact that in single molecules experimentalists were able to measure 
electronic conduction simultaneously with the photon emission^ or Raman response^ is 
very encouraging for the present proposal. What one should monitor is the current / at 
T < Tk with an applied small dc source-drain voltage and subject to a monochromatic FIR 
radiation with tunable frequency. To anticipate, most important for this experiment is that 
the absorption intensities need not be measured. The manner to conduct the experiment 
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and to deduce the relevant parameters can easily be understood by inspecting Fig. [61 

Let us assume that the gate potential Vg is increased, starting from a sufficiently negative 
value (Vg < Vg^i), where the dot level is empty, and / = 0. As soon as the Kondo plateau 
is reached {Vg > Vg^i, Ed < £f), which is signaled by the onset of a current flow / 7^ 0, FIR 
absorption becomes possible by appropriately tuning the frequency u of incoming photons, 
u = uJs{Vg). For ascertaining the resonance, it is not necessary to detect a nonvanishing 
absorption intensity: the resonance will be signaled by the current drop (/ ~ 0), which 
should be observable in a time-resolved dc-transport measurement, because, by absorbing a 
photon, the unpaired electron of the dot will be displaced into electrodes, and the prerequisite 
for the Kondo effect will disappear. In this region, the second signal at the higher frequency 
could hardly be detected, because of its small intensity (see Fig.[T)D), and this is visualized 
by the dashed line in Fig. [61 By further increasing Vg, it will acquire sufficient intensity; 
the representation in Fig. [61 switches from a dashed to a solid line. There, a current drop 
is observed by tuning the photon frequency both to u = uj^iVg) and to u = u!4{Vg). The 
frequencies and U4, which vary linearly with Vg, become closer and closer, and the two 
absorption signals tend to coalesce. Their overlap is perfect {u^ = UJ4) at Vg = V* at the 
point of particle- hole symmetry. Beyond this point, the trend reverses: the lower frequency 
UJ3 decreases while the upper frequency UJ4 increases, and the latter signal eventually becomes 
too weak to produce a current drop (the line switches from solid to dashed). 

Importantly, U and a can be determined from the curves uj^^^^Vg). The former can be 
deduced from the location of the two overlapping absorption peaks in the middle of the 
Kondo plateau, U = 2uj3^4(y*). If there were uncertainties to exactly locate the point of 
the perfect overlapping, one could alternatively use the intersection of the extrapolated u^- 
and u;4-lines. This is a direct determination of U, and not an indirect one, as in low-bias 
dc-measurements, for which the conversion factor a is needed. Moreover, even a can be 
directly obtained from the slope of the curves, a = ^duj3^4/dVg, or, alternatively, by using 
the value of U [a = U/iVg^i — V^,«)]. So, one can even perform a self-consistency test. 
Once an accurate a-value is available, one can use the extension 5Vg of the Kondo plateau 
edges to obtain the parameter F ~ f^/D, which characterizes the finite level width induced 
by the QD-electrode coupling. This is also important, because can also be controlled 
experimentally by varying the gate potentials that form the constrictions.- 
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V, v' V Vo 

FIG. 6: (Color online) Schematic representation of the change in the FIR absorption frequencies 
103^4 and dc current / by varying the gate potential Vg, which permits to deduce the parameters U, 
a, and F. See the main text for details. 

IX. CONCLUSION 

The smaller the size of a quantum dot, the larger is its charging energy. Larger QDs 
possess smaller charging energies (for example, U ~ 64fieV-^), and therefore they could be 
investigated by rf or microwave techniques. However, smaller QDs, as those often used in a 
SET setup, are characterized by considerably larger charging energies (e. g., U ~ 1.9 meV,- 
OT U ~ 7 — SmeV^^), and, consequently, for them the aforementioned techniques cannot be 
directly employed. In the present paper, we have presented theoretical results demonstrating 
that FIR experiments on such SETs, which are feasible, permit to accurately determine the 
charging energy and other important parameters in a direct way. Concerning the FIR 
absorption, three aspects are worth to be mentioned. 

First, we emphasize that, in comparison with other methods, the FIR absorption possesses 
important advantages. It is not affected by parasitic currents due to unavoidable capacitive 
couplings, as it is the case of rf or microwave techniques. Likewise, it is much less challenging 
than photoionization studied recently:'^ in the FIR experiments discussed in the present 
paper one simply needs to determine absorption energies of the order of a few meV with a 
reasonable accuracy, while photoionization requires the determination of ionization energies 
(of the order of the work functions, typically ~ 1 eV) with an accuracy ~ 1 meV.— 

Second, we note that the investigation with the aid of FIR radiation is by no means limited 
to SETs. In nanodevices based on double (or other assembled) QDs, FIR absorption can 
also be used to deduce other relevant parameters,-^ like the interdot electrostatic coupling 
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or V^-Hubbard strength), which are related to important properties of nanostructures (see, 



e.g., Refs. 



23 



and 



25l ). and which cannot be straightforwardly deduced from zero-bias dc- 
conductance data. 

Third, one should emphasize the cross-fertilization between NRG and exact numerical 
diagonalization. Based on a few significant many-body configurations, the latter method 
is very intuitive and allowed us to give a physical content to the NRG numerical findings 
unraveled so far. Conversely, the agreement between the NRG results, valid for infinite 
electrodes, and the exact diagonalization, which can be carried out only for short electrodes, 
demonstrates that the latter is able to make certain valuable predictions that are not affected 
by finite-size effects, as already noted.— >^ 

In addition to the FIR absorption, in the present paper we presented results on the SET 
microwave/rf absorption, which, although preliminary, are interesting in the context of the 
recent experiments revealing the splitting of the Kondo resonance.— We hope to return 
soon to this important issue, which deserves further work. 
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